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0^ : 1 Introduction 

w^ ■ In studying the equation of state (EOS) of ordinary quark matter, the cruial point is 

to treat quark confinement in a proper way. Except the conventional bag mechanism 
(where quarks are asymptotically free within a large bag), an alternative way to 
^ I obtain confinement is based on the density dependence of quark masses, then the 

^ ■ proper variation of quark masses with density would mimic the strong interaction 

between quarks, which is the basic idea of the quark mass-density-dependent model. 
^ , Originally, the interaction part of the quark masses was assumed to be inversely 

C^- i proportional to the density (Fowler et al. 1981; Chakrabarty 1991; Chakrabarty 1993; 

Chakrabarty 1996), and this linear scaling has been extensively applied to study the 
Tj- ■ properties of strange quark matter (SQM). However, this class of scaling is often 

t~^ ■ criticized for its absence of a convincing derivation (Peng 2000). Then a cubic scaling 

^ , was derived based on the in-medium chiral condensates and linear confinement (Peng 

2000). and has been widely used afterwards (Lugones & Horvath 2003; Zheng et 
al. 2004; Peng et al. 2006; Wen et al. 2007; Peng et al. 2008). But this deriving 
procedure is still not well justified since it took only the first order approximation of 
the chiral condensates in medium. Incorporating of higher orders of the approximation 
would nontrivially complicate the quark mass formulas (Peng 2009). In fact, there 
are also other mass scalings in the literatures (Dey et al. 1998; Wang 2000; Zhang et 
al. 2001; Zhang & Su 2002; Zhang & Su 2003). 

Despite the big uncertainty of the quark mass formulas, this model, after all, is no 
doubt only a crude approximation to QCD. For example, the model may not account 
for quark system where realistic quark vector interaction is non-ignorable. However, 
we can not get a general idea of how the strong interaction acts from the fundamental 
theory of strong interactions in hand, i.e. QCD. Until this stimulating controversy 
is solved, we feel safe to take the pragmatic point of view of using the model. This 
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work does not claim to answer how Nature works. However, it may shed some hght 
on what may happen in interesting physical situations. In this respect, the quark 
mass-density-dependent model has been, and still is, an interesting laboratory. 

The aim of the present paper then, is to study in what extent this scaling model 
is allowed to study the properties of SQM. To this end, we treat the quark mass 
scaling as a free parameter, to investigate the stability of SQM and the variation 
of the predicted properties of the corresponding strange stars (SSs) within a wide 
scaling range. Furthermore, we try to demonstrate the general features of SSs related 
to astrophysics observations, whatever the value of the free parameters. 

The paper is organized as follows. In Section 2 we describe the formalism applied 
in calculating the EOS of the SQM in the quark mass-density-dependent model. In 
Section 3 we present the structure of the stars made of this matter, including mass- 
radius relation, spin frequency, electric properties of the quark surface. Finally in 
Section 4 we address our main conclusions. 

2 The Model 

As usually done, we consider SQM as a mixture of interacting -u, rf, s quarks, and 
electrons, where the mass of the quarks niq {q = u, d, s) is parametrized with the 
baryon number density rib as follows: 



rrig = irigo + mi = m^o + — , (1, 
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where C is a parameter to be determined by stability arguments. The density- 
dependent mass niq includes two parts: one is the original mass or current mass 
mqo, the other is the interacting part mi. The exponent of density x, i.e. the quark 
mass scaling, is treated as a free parameter in this paper. 

Denoting the Fermi momentum in the phase space by z/j {i = u,d,s,e~), the 
particle number densities can then be expressed as 






-.-.. 7.^-^/ PMP-S, (2) 



and the corresponding energy density as 



^ = E^^"vW^P'dp. (3) 

The relevant chemical potentials /i„, Hd, Us, and He satisfy the weak-equilibrium 
condition (we assume that neutrinos leave the system freely): 

l^u + l^e=l^d, Hd = f^s- (4) 



For the quark flavor i we have 
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We see clearly from Equ. ([5]) that since the quark masses are density dependent, the 
derivatives generate an additional term with respect to the free Fermi gas model. 
For electrons, we have 
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The pressure is then given by 

P = -e + Y^ fiiUi 



=u,d,s \ J / 
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with VLq being the free-particle contribution: 
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The baryon number density and the charge density can be given as: 



nh = -{riu + nd + ns), 
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Qq = 7T"-« - 7;TT-d - -TT-s - 'T-e- 



(10) 
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The charge- neutrality condition requires Qq = 0. 

Solving Equs. @, f lTO|) . flTTl) . we can determine Uu, Ud, Ug, and rig for a given total 
baryon number density rib- The other quantities are obtained straightforwardly. 



In the present model, the parameters are: the electron mass nie = 0.511 MeV, 
the quark current masses rriuo, fnao, rrisQ, the confinement parameter C and the quark 
mass scaling x. Although the light-quark masses are not without controversy and 
remain under active investigations, they are anyway very small, and so we simply 
take niuo = 5 MeV, rrido = 10 MeV. The current mass of strange quarks is 95 ± 25 
MeV according to the latest version of the Particle Data Group [19] 
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Figure 1: The stability window of the SQM at zero pressure with the quark mass 
scaling parameter x = 1/3,1,2. The stability region (shadow region), is where the 
energy per particle is lower than 930 MeV and two-flavor quark matter is unstable. 

We now need to establish the conditions under which the SQM is the true strong 
interaction ground state. That is, we must require, at P = 0, E/A < M(^^Fe)c^/56 = 
930 MeV for the SQM and E/A > 930 MeV for two-flavor quark matter (where 
M(^^Fe) is the mass of ^^Fe) in order not to contradict standard nuclear physics. The 
EOS will describe stable SQM only for a set of values of (C, mso) satisfying these two 
conditions, which is given in Fig. 1 as the "stability window" . Only if the (C, rriso) 
pair is in the shadow region, SQM can be absolutely stable, therefore the range of C 
values is very narrow for a chosen rrisQ value. 

We then illustrate in Fig. 2 the density dependence of mi with the quark mass 
scaling x = 1/10, 1/3, 1, 3. The calculation is done with rriso = 95 MeV and C values 
corresponding to the upper boundaries deflned in Fig. 1 (the same hereafter), that 
is, the system always lies in the same binding state (for each x), i.e, E/A = 930 
MeV. We presented those C values in the last row of the Table. 1. Clearly the quark 
mass varies in a very large range from very high density region (asymptotic freedom 
regime) to lower densities, where conflnement (hadrons formation) takes place. It is 
compared with Dey et al.'s scaling (dash-dotted) [6]. 

3 Results and Discussion 



The resulting EOSs for the SQM are shown in the left panel of Fig. 3 with x = 

1/10, 1/5, 1/3, 1, 2, 3. Because the sound velocity v =\ dP/dp \ ^/^ should be smaller 




Figure 2: The density dependence of mi with the quark mass scahng parameter x = 
1/10, 1/3, 1, 3. The calculation is done with tjIso = 95 MeV and C values presented 
in the last row of the Table. 1 (see text for details). It is compared with Dey et al.'s 
scaling (dash-dotted line) (Dey et al. 1998). 



than c (velocity of light), unphysical region excluded by this condition has been 
displayed with scattered dots. For the x values chosen here, they have quite different 
behavior at low density, basically falling into two sequences. At small scalings {x = 
1/10, 1/5, 1/3) the pressure increases rather slowly with density; while the curve turns 
to rapidly increase with density at relatively large x values (x = 1,2,3). They cross 
at e ~ 800 MeV fm ~^, then tend to be asymptotically linear relations at higher 
densities, and a larger x value leads to a stiffer EOS. 

This behavior of EOSs would be mirrored at the prediction of mass-radius relations 
of the corresponding SSs, as is shown in the right panel of Fig. 3. For the first 
sequence, the maximum mass occurs at a low central density (as shown in Table. [1]), 
so a higher maximum mass is obtained due to a stiffer EOS, and with the increase of 
X value, the maximum mass is reduced from 1.78M0 at x = 1/10 down to 1.61 Mq 
at X = 1/3; While we observe a slight increase of the maximum mass with x value 
for the second sequence: from 1.56M0 at x = 1 up to 1.62 Mq at x = 3. Anyway 
the resulting maximum mass lies between I.5M0 and I.8M0 for a rather wide range 
of X value chosen here (0.1 - 3), which may be a pleasing feature of this model: well- 
controlled. To see the region of stellar parameters allowed by this model, we plot in 
Fig. 3 also the M(R) curves for lower boundaries defined in Fig. 1 with x = 1/5, 1/3, 1 
(grey lines in the right panel). 

Moreover, the radii invariably decrease with x value. In addition, we employ the 
empirical formula connecting the maximum rotation frequency with the maximum 
mass and radius of the static configuration [S], and present also the maximum ro- 
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larger x value results in a larger maximum spin frequency, SSs with a; = 3 can rotate 
at a frequency of 2194 rad s~^. More detailed results can be found in Table. [TJ 
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Figure 3: The SQM EOS and mass-radius relation of SSs with x = 
1/10,1/5,1/3,1,2,3. M(R) curves for lower boundaries defined in Fig. 1 with 
X = 1/5,1/3,1 are also presented (grey lines in the right panel). Contours of the 
maximum rotation frequencies are given by the light grey curves (Gourgoulhon et al. 
1999). 



In addition, the surface electric field should be very strong near the bare quark 
surface of a strange star because of the mass difference of the strange quark and the 
up (or down) quark, which could play an important role in producing the thermal 
emission of bare strange stars by the Usov mechanism (Usov 1998; Usov 2001). More- 
over, the strong electric field plays an essential role in forming a possible crust around 
a strange star, which has been investigated extensively by many authors (for a re- 
cent development, see Zdunik et al. 2001). Also it should be noted that this electric 
field has some important implications on pulsar radio emission mechanisms (Xu et al. 
2001). Therefore it is very worthwhile to explore how the mass scaling influences the 
surface electric field of the stars, and possible related astronomical observations in 
turn may drop a hint on what the proper mass scaling would be. Adopting a simple 
Thomas- Fermi model, one gets the Poisson's equation (Alcock et al. 1986): 



^_/ tiv'-v') ^<0' 



dz^ 



^y3 
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z > 0, 



(12) 



where z is the height above the quark surface, a is the fine-structure constant, and 
V^ /{^Ti'^h c^) is the quark charge density inside the quark surface. Together with the 
physical boundary conditions {z — )■ — oo : V ^ Vq, dV/dz — )■ 0; ^ — )■ +oo : V — )■ 
0,dV/dz — i- 0}, and the continuity oi V at z = requires V{z = 0) = 3Ki/4, the 
solution for z > finally leads to 



V 



3K 



'-V^z + 4 



(for z > 0). 



(13) 



The electron charge density can be calculated as V /{Sn h c ), therefore the number 



density of the electrons is 



n. 
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and the electric field above the quark surface is finally 



E 




^v^-z + Ay 



(15) 



which is directed outward. 

It is shown in Fig. 4 (take a; = 1/3 for example) that although the electric field near 
the surface is about 10^^ V cm~^, the outward electric field decreases very rapidly 
above the quark surface, and at z ~ 10~^ cm, the field gets down to ~ 5 x 10^^ 
V cm~^, which is of the order of the rotation-induced electric field for a typical 
Goldreich- Julian magnetosphere. To alter the mass scaling mainly has two effects: 
First, it affects a lot the surface electric field, and a small scaling parameter leads to 
an enhanced electric field. The change of electric field would be almost a order of 
magnitude large (from 10^^ V cm~^ to 10^^ V cm~^), which may have some effect on 
astronomical observations. Second, a larger scaling would slow the decrease of the 
electric field above the quark surface. 




Figure 4: The electric field above the quark surface with x = 1/3, 1, 2, 3. 



4 Conclusions 

In this paper, we investigate the stability of SQM within a wide scaling range, i.e. 
from 0.1 to 3. We study also the properties of the SSs made of the matter. The 
calculation shows that the resulting maximum mass always lies between 1.5Mq and 
I.SMq for all the mass scalings chosen here. Strange star sequences with a linear 
scaling would support less gravitational mass, a change (increase or decrease) of the 
scaling parameter around the linear scaling would result in a higher maximum mass. 
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Radii invariably decrease with the mass scahng; and then the larger the scaling, the 
faster the star rotates. In addition, the variation of the scaling may cause an order 
of magnitude change of the surface electric field, which may have some effect on 
astronomical observations. 
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M/Mq 


1.78 


1.61 


1.58 


1.56 


1.61 


1.62 


R/km 


13.2 


9.38 


8.75 


8.10 


7.97 


7.89 


ric/riQ 


4.35 


7.88 


8.88 


10.1 


10.2 


10.3 


Qma^/rad s~^ 


1066 


1691 


1860 


2072 


2159 


2194 


C/MeV fm"^^ 


199.1 


126.8 


104.1 


69.5 


41.7 


28.8 



Table 1: Calculated results for the gravitational masses, radii, central baryon densities 
(normalized to the saturation density of nuclear matter, uq = 0.17 fm~^), and the 
maximum rotational frequencies for the maximum-mass stars of each strange star 
sequence. The calculation is done with niso = 95 MeV and C values presented in the 
last row of this table. 
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